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We develop a formalism for constructing the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) ma- 
trix and neutrino masses using an expansion that originates when a sequence of heavy right handed 
neutrinos are integrated out, assuming a seesaw mechanism for the origin of neutrino masses. The 



o 

fS| expansion establishes relationships between the structure of the PMNS matrix and the mass differ- 

5P ences of neutrinos, and allows symmetry implications for measured deviations from tri-bimaximal 

.^^ form to be studied systematically. Our approach does not depend on choosing the rotation between 

■r^lj- the weak and mass eigenstates of the charged lepton fields to be diagonal. We comment on using this 

expansion to examine the symmetry implications of the recent results from the Daya-Bay collabo- 

r~| ration reporting the discovery of a non zero value for 6113, indicating a deviation from tri-bimaximal 

1^ form, with a significance of 5.2 tr. 

Oh 
cn i. introduction 

> 

o 

^ The standard theory of neutrino oscillations, where three mass eigenstate neutrinos differ from their 

interaction eigenstates leading to the observed neutrino oscillations, is consistent with current experimental 
O data. The amplitude of the oscillations among various neutrino species is related to the misalignment of the 

^ interaction and mass eigenstates of the neutrinos, characterized by the Pontecorvo-Maki-Nakagawa-Sakata 

• ^H (PMNS) matrix |[T]|3 . Taking the relation between the interaction (primed) and mass (unprimed) eigenstates 

;-( to be given by f'^^m = ^(/) L/R) fi/R, where f = {u,e, /j.,- ■ ■ }, the PMNS matrix is given by 

UpMNS=U\e,L)U{v,L). (1) 

The UpMNS matrix can be parameterized in terms of three angles ^12,^13, ^23 and three CP violating phases 
(5, ai^2- Defining Sij = sin Oij and cij = cos Oij with the conventions < Oij < ti 12, and < (5, a\^2 < 27r 
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a general parameterization of this matrix is given by 
/l o\ / 



UPMNS 



C23 S23 
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' Ci2 Si2 
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Cl3 5136 

1 

\ -5136-*'^ Ci3 / 

where P = diag(e*"i'^, e*"^'^, 1) is afunctionof the Majorana phases, present if the right handed neutrinos 
are Majorana, while (5 is a Dirac phase. This latter phase can contribute in principle to neutrino oscillation 
measurements, while the Majorana phases cannot. Recent global fit results on neutrino mass differences and 
measured mixing angles using old/new reactor fluxes are given in Ref. [31, (with new reactor flux results in 
brackets): 
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-{m\+m\)/2 and ATTi^g > 0, (< 0) corresponds 



to a normal (inverted) mass spectrum. This pattern of experimental data is perhaps suggestive of a PMNS 
matrix that has at least an approximate "tri-bimaximal" {TB) form |]41. Fixing sin^(6'i2) = 1/3 and ^13 = 
the TB form is 



UPMNS ~ T^TB 
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(4) 



for a particular phase convention. 

The structure of this matrix could be fixed by underlying symmetries. In attempting to determine such 
an origin of this matrix, the "flavour" basis where one assumes U{e, L) = diag(l, 1, 1) is frequently used. 
When this assumption is employed the relationship between the weak and mass neutrino eigenstates is 
identified with the Upmns^ i-C- v[ = (UpMNs)ij ^j for hj = 1)2,3. There have been many attempts 
to link the approximate TB form of the neutrino mass matrix to symmetries of the right handed neutrino 
interactions in this basis, see Ref. ||5] El for a review. Recent experimental results provide evidence for 
deviations from this TB form. Evidence for sin^(6'i3) 7^ in global fits is reported to be > 3(7 in Ref. ||3l 
at this time. As this paper was approaching completion, the discovery of non- vanishing 6*13 was announced 



by the Daya Bay Collaboration Q with a reported value of 

sin2(2 (9i3) = 0.092 ± 0.016(stat) ± 0.005(syst), (5) 

corresponding to 5.2a evidence for non zero ^13. It is reasonable to expect further speculation about the 
origin of the deviation from TB form in light of this result, where again the flavour basis will be frequently 
assumed. 

There is no clear experimental support for assuming that U{e,L) = diag(l, 1, 1). This choice can 
be motivated by an ansatz related to the origin of the approximately diagonal structure of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, and a further ansatz on the relation between U'^{d, L) and U'^{e, L), 
see Ref. ||5l for a coherent discussion on this approach. This choice can also be justified with model building 
in principle, see Ref. |[8l for an example. Conversely, in grand unified models frequently associated with the 
high scale involved in the seesaw mechanism, the quark and lepton mass matricies can be related in such 
a manner that the flavour basis cannot be chosen, or at least, the choice of the flavour basis can be highly 
artificial. 

It is of interest to have a formalism for neutrino phenomenology that is as robust and basis independent 
as possible. Clearly linking symmetries to the form of the PMNS matrix requires that the physical conse- 
quences of a symmetry are not dependent on a basis that can be arbitrarily choosen for Z^(e, L), such as the 
flavour basis. In this paper, we develop a perturbative approach to the structure of the PMNS matrix as an 
alternative to symmetry studies that are basis dependent. ^ Using this approach we show a basis independent 
relationship between the eigenvectors that give U{v, L) and U{e, L) corresponding to TB form at leading 
order in our expansion. We also examine the prospects of relating patterns in the data in low scale measure- 
ments of neutrino properties with high scale flavour symmetry breaking, illustrating how our approach can 
be used to study the recent reported discovery of non zero ^13 reported in Ref. Q. 

II. CONSTRUCTING A FLAVOUR SPACE EXPANSION FOR NEUTRINO PHENOMENOLOGY 

In this section we develop a formalism linking the measured differences in the neutrino mass eigenstates 
with the structure of the PMNS matrix, which we will refer to as a flavour space expansion (FSE). There 
is an experimental ambiguity in the measured mass hierarchies at this time. The neutrino mass spectrum 
can be a normal hierarchy (7713 > m2 ^ mi) or an inverted hierarchy {m^ < mi ^ 7712). In our initial 
discussion we will assume a normal hierarchy. The formalism can be reinterpreted for an inverted hierarchy. 



' This is a modem implementation in the neutrino sector of an old idea of relating the mixing matricies of the standard model 
(SM) to the measured quark or neutrino masses, for pioneering studies with this aim see Ref. l9lll4tfT6l . 



A. Review of Seesaw the Mechanism 

Recall the standard seesaw scenario IITtII . with three right handed neutrinos Nm? We use the notation 
i = {vl, cl) for the left handed SU(2) doublet field and cr for the SU(2) singlet lepton field carrying 
hypercharge. These fields carry flavour indicies i,j? We also define H = iT2 H* , where H is the Higgs 
doublet of the SM, with (-ff^) = (0, u/\/2). Then the lepton sector of the Lagrangian in the seesaw scenario 
is the following 

C = i N'rJN'j,, -^NgMl, N'^, - 4. H {y'E\, e'^^ - iVJj, {y',\, H^ l'^, + h.c. , (6) 

here we have defined Nf^^ = CN'^^ with C the charge conjugation matrix. There is freedom to rotate to the 
mass basis by introducing the unitary rotation matricies U{f, L/R), so that 

N'jii = U{N, R)^j Nrj, 4, = U{L, L)i, hj, e'ffi = U{e, R)ij cr,. (7) 

The kinetic terms are unchanged by these rotations, and one is free to further rotate to a basis in the flavour 
space of M, tje, Vv defined by the set of all possible M, i/e, Vv related through 

M' -^ U{N,RY MU{N,R)\ 

v'e ^ U{L,L)yEU{e,R)\ (8) 

y', ^ U{N,R)y,U{L,L)^. 

We choose to make the initial rotation to a basis where M is diagonal real and nonnegative. This fixes 
U{N,R). Initially the Majorana mass matrix has three complex eigenvalues Ma^b,c = f]a,b,c |-^a,6,c| = 
''la,b,c {fJ'a,b,c), here 7]a^b,c are the Majorana phases. Working in this basis llT8l420l shifts any Majorana phases 
into the y,^ matrix and Ni = ^/rjiNRi + \/Vi ^m- Integrating out the heavy Ni one obtains the dimension 
five operator ||2T1 

C^ = \{HUi)Cij{HUj) + h.c., (9) 

with the matrix Wilson coefficient Cij = (yj^ rjM~^ yv)ij- The key observation that we use in this work is 
that when integrating out the Ni in sequence"* a perturbative expansion for neutrino phenomenology that is 

^ Our initial discussion will largely follow Ref. 11814201 . 

^ These flavour indices will run over a, h, c for the three Ni. We use this notation to distinguish these flavours from the measured 

mass differences of the physical eigenstates (7711,2,3) and also the index on the Yukawa vectors of each N[, which run over 

1,2,3. 
* We distinguish in this paper right handed neutrinos that are in a mass diagonal basis, reserving the notation Ni for such states, 

compared to NRi states, which are not mass diagonal in general. 



related to a hierarchy in the magnitude of the contributions to Cij can be constructed. Note that this is also 
the key point in the sequential dominance idea of Refs. ifTOUTBl . however, the formalism we will develop is 
distinct from these past results. 

B. Perturbing the Seesaw 

1. Model Dependence of Perturbing the Seesaw 

The FSE we develop depends on the mass spectrum and the Yukawa couplings of the Ni. Also, in 
principle, the flavour orientation^ effects the quality of the FSE. Experimentally, at this time, all of these 
Lagrangian parameters are individually unknown, thus we will be forced to assume that a perturbative 
expansion of the form we employ exists. Our approach is to view the seesaw Lagrangian given in Eq. 
(6) as an effective theory, and we do not attempt to uniquely identify and restrict ourselves to a particular 
UV physics that dictates the low energy parameters in this effective theory in this paper. However, the 
formalism we develop allows perturbative investigations of the neutrino mass spectrum and mixing angles 
in many scenarios and is in fact quite general. It would be surprising if all of the unknown parameters 
conspired to forbid an expansion of this form from being present. We seek to illustrate this point in this 
section, by demonstrating a simple scenario where the FSE would be of some utility. 

There are of course cases when the formalism we outline cannot be used. When the ^ia,b,c and Yukawa 
couplings of the Ni ai^e both highly degenerate a FSE cannot be used.^ This is also the case if a non- 
degenerate pattern of the ^j.a,h,c and the Yukawa couplings are such that the perturbations to the neutrino 
masses are similar in magnitude. However, as the iia,b.c depend on the matching of the effective theory to 
UV physics, while the Yukawa couplings are dimension four and are not UV sensitive, this would require 
tuning the physics of different energy scales. Unless model building is used to relate the masses and Yukawa 
couplings of the Ni so that the contributions to the vl masses are the same, a FSE is expected quite generally, 
and can be related to neutrino mass differences, as we will show. 

One can study the expected spectrum of left handed masses due to Eqn. (|6]) in generic scenarios. Without 
assuming other interactions beyond the SM, the N'^^ are not distinguished by any quantum number. As a 
result, the non-diagonal mass matrix that is the coefficient of the operator matrix Oij = N^,-^ N'^- given 
by Mj' is naively expected to have entries that are all similar in magnitude. As the operator is dimension 



^ By the flavour orientation of the Yukawa coupUngs, we mean the orientation of the Yukawa coupling vectors x, y, z with respect 

to the leading order eigenvectors pa,b,c', see the next section for definitions and further details. 
^ An exception to this statement is when the flavour orientation is such that an expansion is present due to the geometry in flavour 



space, we discuss an example of this form in Section II D 



three, the mass matrix is expected to be proportional to the highest scale UV physics (violating lepton 
number) that was integrated out leading to this effective theory. We denote this scale by Mq and the naive 
non-diagonahzed mass matrix in this case as 



/ 



M'- ~ 



Mo 



1 1 1 



\ 



(10) 



1 1 1 
\l 1 1/ 

with eigenvalues //? = {3 Mq, 0,0}. The vanishing of the two eigenvalues of M[- can be lifted by interac- 
tions of the N'^^. These interactions can be dictated by beyond the SM quantum numbers assigned in UV 



model building. If the mass matrix is still approximately degenerate, as in Eq.( 10 1, then a hierarchy of the 
Ni masses is still expected. An example of a small breaking where other interactions do not need to be 
assumed is given by the orientations of the N'^^^ in flavour space. Rotating to the lepton diagonal mass basis, 
these interactions give loop corrections 
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suyui^ ^ (^^)y_^p-^-^ log 
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(11) 



Logarithmic corrections of this form are required to cancel the renormalization scale dependence of the pole 
masses of the N'^^. Including such corrections leads to M^' ~ Mq (1 + eij). These corrections split the 
mass spectrum; diagonalizing one finds 



^1 = Mo 3 + 
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ij ^ij 



fJ-2 



Mn 
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Kj 



fi3 = MoO{e^). 



(12) 



A hierarchical spectrum of ^j is expected with the pattern (/ii, ;U2, /^s) ~ Mo(l, e, e^). 

In the case of degenerate /Uj due to the M^- mass matrix not conforming to these generic expectations, 
the expansion can follow from a hierarchy in the Yukawa couplings of the Ni - such as the hierarchical 
pattern of Yukawa couplings in the SM. Lastly, the suitability of the expansion can follow only from the 
orientation in flavour space of the Yukawa coupling vectors of the Ni. In summary, the appropriateness of 
the FSE is clearly model dependent, but we expect it to be broadly applicable in realistic models. 

C. Developing the Seesaw Perturbations 

Consider supplementing the SM field content by a single right handed neutrino' Na and an interaction 
term that couples it into a linear combination of £^. The coupling is fixed by the complex Yukawa vector 



^ The subscript in A'^; will run over labels a, b, c; for now we are just considering the single spinor field A^a- 



x^ = {xi,X2,x^} in flavour space whose form can be constrained by a flavour symmetry, but is here left 
arbitrary. We absorb the overall Majorana phase into this vector. The relevant Lagrangian terms^ are given 
by 

Ca = -^NaNa-XpNaHUl + h.C. (13) 

The right handed neutrino can be integrated out, giving for the left handed neutrino mass matrix a nonzero 
eigenvalue 

M = - — xx^ = U{u, LY diag(0, 0, m^) U{v, L)^ . (14) 

The matrix U{u,L) is the matrix {pY , pi* , pi*) of normalized (column) vectors p*f^^ that solve 
Aip* = mp, with m real and non-negative.^ These vectors are also eigenvectors of A^^ 7W = 
{x*x^) {xx'^) = {x^x)x^x'^ with eigenvalues m?; their complex conjugates, p are eigenvectors of 
TWTW^ = (xx-^) (x^af^) = {x^x)xx^, also with eigenvalues m"^. One finds the leading eigenvector and 
eigenvalue 

Pa = x/\x\, ma = V^{x^x)/2pa- (15) 

Here pa is real and non-negative due to the initial flavour basis choice that fixed U{N, R), while U{v, L) 
and x are in general complex.^" The remaining two eigenvectors pj, c are such that {pa\pb) = Pa Pb = 0, 
{Pa\pc) = 0. These eigenvectors will lead to the remaining (smaller) mass eigenvalues, and will perturb 
the leading eigenvalue and eigenvector and thus the U{v,L) matrix. This is the expansion we seek to 
exploit. Consider the perturbation that generates the second eigenvalue of the neutrino mass matrix due to 
a second right handed neutrino N^, which we define to couple into a linear combination of the i^ given 
by 2/^ = {yi, 2/2, 2/3}- One obtains a second eigenvalue in the left handed neutrino mass matrix so long as 

y [f X. The perturbation of M. M!^ is given by 

44 4 

{M + 5M){M + 6MY = ^xiF{x*x^) + -^^{xx^){y*y^) + -^^{yf){x*x^), 

^Pi 4:Pa Pb ^Pa Pb 

,,4 



+ ^yfiy'y')- (16) 

^Pb 

At leading order the pb^c have degenerate (vanishing) eigenvalues. One is free to rotate to a chosen basis in 
these vectors. We rotate to a basis in these vectors such that the following perturbation vanishes for pc 

{Pc\{yy'^){y'^y*)\Pc)=^. (17) 



^ We have chosen U{N, R) as discussed in Sectionplto rotate to a mass diagonal basis for Ni and absorbed this rotation into a 
redefinition of the Yukawa matrices y^. The Yukawa vectors x, yand i* introduced below are given in this basis. 

' Here we use the notation that the eigenvectors run over the flavour indicies a, fe, c to denote that they are associated with inte- 
grating out each of the Na.t.c at leading order in the corresponding mass eigenvalues. 
'" This choice is allowed by the unitary flavour transformations that are symmetries of the kinetic terms, Eqs. (6-8). 



With this choice pc retains a vanishing eigenvalue when Ni, is integrated out. The eigenvector pc should be 
orthogonal to pa oq x. A normalized eigenvector basis at leading order is then 

_^ X* X (y X x) ^ y* X X* 

Pb = 1^1 1 ^ — =r' P'^ = "R — =r- (1^) 

|x| \x X y\ \x X y\ 

Now we can determine the perturbation on the leading order eigenvectors and eigenvalues when Nf, is 
integrated out. The corrections to the eigenvalues and eigenvectors using perturbation theory are given by 

^^ {pi\M6M^ + 5M M^Pi) ^ 

°Pj = = Z^ 2 2 Pi^ 

Smj = {pi\M6M^ + 6MM^ + 6M6M^\pi). (19) 

Nonzero eigenvalues are obtained for 7715 c at second order in the expansion due to the orthogonal basis 
vectors causing the leading perturbations to each of these masses to vanish. The leading perturbation to the 
eigenvectors and ma is first order in the expansion. We retain the leading perturbation on the eigenvectors 
and the leading and subleading effects on the masses to obtain nonzero eigenvalues. We find the following 
for the perturbations 
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c- 2 {y\pb){x-y )\x\ ^ 

^Pb = l^ab T^ Pa, 

5pc = 0, 6m'i = 0. (20) 

Here (5m?- = m,f — m|, pf, = v'^/{4,pipj) and cosOxy = \x* ■ y\/\x\\y\ is a measure of (the cosine 
of) the angle between the x, y Yukawa vectors. The masses are evaluated to the appropriate order in 
the perturbative expansion and the eigenvector perturbations are in general complex. Note that for the 
phenomenology of the Upmns matrix that we will study it will be sufficient to only retain the leading 
perturbation, while when studying the mass spectrum the leading and sub-leading perturbations should be 
retained. 

Finally integrate out the third right handed neutrino Nc with Majorana mass pc, which couples into a 
linear combination of the £^ dictated by zT = {zi,Z2,z^}. The resulting eigenvector perturbations are 



2 {Pb\z){z-X*)\x\ ^ 2 {Pc\z){z-X *)\x\ 

'"^ 5ml, ''' ^"^ 5mh 



52Pa = Pac ^^ Pb + fJ-ac ^^ Pc, 

^ab 



2 {z\Pb)ix-Z *)\x\ 

''"' 5ml 



(>2Pb = ^J■ac ^^ Pa, 

\a 

e ^ 2 {z\Pc){x ■ Z*)\x\ ^ 

02pc = f^ac r^^ Pa- (21) 



5m1„ 



The mass perturbations are 

52ml = 2 fil Re [{z\x) {x ■ z')] + fil \z\^ cos^ 9, 



XZl 

2 _ ,,2 l/;T*l;rM2 I --12 



hm^, = ^icc\{pb\^\ \z\ , 

52ml = ^il,\{pc\^\''\z\\ (22) 

The measured masses of the neutrino's are related to these perturbations as 

m\ = ml + 5ml + 52ml 

ml = 5ml + 52ml 

mi(j = 52m,^. (23) 

It is interesting to note that a normal hierarchy emerges quite naturally from the FSE as the leading neutrino 
mass rua receives corrections at linear order to its mass, while the remaining masses only receive corrections 
at second order in the perturbations. 

It is also important to note that this formalism does not require a hierarchy of the form 52m!f <^ 5ml, 
only 52m,l, 5m,l <^ ml is required. Expanding on this important point in more detail, it is not required that 
the perturbation due to integrating out Nf, is larger than the perturbation due to integrating out Nc- Only that 
the effect of integrating out each of these right handed neutrinos perturbs the initial mass matrix — which 
is dominated by integrating out the initial right handed neutrino Na- The existence of these perturbations 
are not necessarily direct statements on the relative size of the /ij as we discuss in more detail in the next 
section. 

D. Inverted and Normal Hierarchy and Flavour Space 

For a normal hierarchy (with notation 771,3 > "t^2 ^ m.i) we identify (1,2,3) = {C,B,A) in the 
equations above. The difference between the normal and inverted (7713 < ttii ^ 7712) case appears in the 
relative size of the 5mi and 52m,i and the identification (C, B, A) = (3, 1, 2) for an inverted hierarchy. An 
inverted hierarchy requires non generic perturbations in the FSE, or one can trivially modify the expansion 
to only perturb when Nc is integrated out. The size of drrii and 52m,i depends on the hierarchy in iJ.a,b,c, the 
magnitude of the Yukawa vectors, and also the orientation in flavour space of the vectors x, y, z. In this 
section, we will discuss the size of the perturbations of the FSE in light of neutrino mass data in a scenario 
where the perturbations are dictated primarily by a hierarchy in fj,a,b,c- We will then discuss the case where 
the expansion arises primarily due to the orientation of the Yukawa coupling vectors in flavour space. 



10 



1. The expansion without flavour alignment 

We can examine the quality of the expansion by comparing to the measured mass differences in neutri- 
nos. Consider the case that the size of the perturbations is generic in the sense that the Yukawa coupling 
vectors are taken to be 0(1) with the orientation in flavour space not significantly effecting the quality of 
the expansion. In this case, the expansion follows from the relative size of the /ij, i.e. ^a < l^b < Re- 
consider the generic case in a normal hierarchy. Using the FSE and retaining the dominant term 

Amis = rnl + 6ml + 62ml - {6ml + 62ml) ^ ml, (24) 

while in the same manner 6m,l r^ Attt-I^ and the expansion requires 62m,'^ < Am32, Attt-I^. Due to this, the 
expansion requires ^^121^/(2 fic) < v^"^ ~ 0-05 eV. This condition is consistent with current bounds 
on the absolute neutrino mass scale, with a 95% C.L. bound of X]m,i^ = 0.28 eV quoted in Ref. ll22]| . 
assuming ACDM cosmology. It is also consistent with current bounds from Tritium /3 decay experiments 
Il23l which quote m{iye) < 2 eV at 95% C.L. 

Expressing this condition in terms of the high mass scale of the Nc integrated out, fJ-c/\z\'^ ~ 10^^ GeV. 
Generically one expects the mass scale of the right handed neutrino operator to be the largest scale integrated 
out that violated L number, and this condition for the lightest neutrino is clearly consistent with naive 
expectations of Mq ~ Mpi. 

2. The expansion with flavour alignment 

Now consider the case where the perturbative expansion follows from the flavour orientation of the 
X, y, z vectors primarily. An example where this is the case is when the threshold matching onto the 
UV physics is such that the Wilson coefficient matrix of Oij yields a mass matrix with nearly degenerate 
eigenvalues. This occurs for example when 

^ l+e e e ^ 

M4. ~Mo e l + e e • (25) 

In this case, a nearly degenerate mass spectrum of the Nji follows, f^a — fJ-b — A^c- Generating a normal 
or inverted hierarchy if one also has \x\ ~ \y\ ~ \z\ requires more precise alignments in flavour space and 
allows a geometric interpretation of the measured neutrino mass spectrum. In the FSE, the tree level masses 
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of the SM neutrinos are 



m\ 



m\ 



mc 



liaa\A +2^a6Fl 1^1 ^os »n,y + 2[i^^x\ \z\ COS 6^^^ + /if,j,|y| cos Br^y^[i^^\z\ cos- 



2|->|4 2/1 I 2|^4 2/1 

/^fefclyl COS Qp^y^i^cAA COS 6lp^, 



/^ccl^^ COs2 6'p^2 



(26) 



Consider the case that all of the /x? ■ are similar in magnitude and | x | > | y | ~ | z| so that the mass spectrum 
is primarily dictated by the orientation of the Yukawa vectors in flavour space. An example of an inverted 
or normal hierarchy is shown in Fig. (1) in this case. 
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FIG. 1: A geometric interpretation of an inverted or normal hierarchy when /x^ ~ /if, ~ /i^ and \x\ ^ \y\ ~ \z\. 
Normal hierarchy on the left, inverted on the right. 



E. UpMNs and Flavour Space 



Now, assuming that a FSE expansion exists, let us consider its utility in examining the form of Upmns- 
The rotation matrix h({v, L) with v* = [pi + Spl + 8291)* is given by U{u, L) = (iJ^, v^, v{). The charged 
lepton mass matrix after electroweak symmetry breaking is given by A^g = vyE/\p2, and diagonalized 
by" 



U{e, R) A^e^ (e, L) = diag{me, m^jirij 



(27) 



" The hierarchy of the charged lepton masses can be used to organize an expansion of W(e, L) in the same manner by constructing 
MlMe in principle. Conversely, in principle, one could employ an ansatz that the hierarchy in these mass eigenvalues could be 
related to an expansion of W(e, R). As such we do not employ a FSE on U{e, L^ . This ambiguity also limits the application 
of a FSE to the CKM matrix. In this manner, the expansion we employ is most useful for expanding W(/^, L) when a seesaw 
mechanism is the origin of the smallness of the neutrino masses. 
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We take Z^(e, L)'^ = ( cfs^, cj2 , cTi^ J , with Bi the orthonormal column eigenvectors diagonalizing AieMl- 
The expanded Up^NS ^^ *^^ "^^ ^^''"^ 



/ 



iY 



PMNS 



+ 0{5pi,52Pi). 



(28) 



173 • 0-3 f 2 • CJ3 Vi ■ 0-3 
^3 • '?2 V2- (72 t^i • (T2 

\v^-ai V2- Ox vi-ai j 

P3- 03 P2- 03 pi • <J3 \ 
= P3 • '?2 P2 • 02 Pl • <?2 

\ P3 • <?! P2 • i?l Pl- Oi J 

This result makes clear that the first two right handed neutrinos that were integrated out in the FSE contribute 
to the leading order structure of the PMNS matrix. The leading order of the pi eigenvectors only depended 
on X, y. As we have discussed, these neutrinos can be integrated out in sequence, and for Yukawa couplings 
that are degenerate the lightest two Ni can lead to the largest mass eigenvalues of the vl. In a mass 
degenerate case, they are more strongly coupled to the i^j. The unitarity of the Upmns matrix is ensured 
when an expansion of this form is employed by normalizing the Vi eigenvectors order by order in the 
expansion. Conversely, if the FSE is used when the expansion parameters are not small, the constructed 
UpMNS matrix is not necessarily unitary. 

It is important to note that, in general, it is the relationship between the eigenvectors that determines 
UpMNS- Exact TB form at leading order in the FSE has a simple interpretation. It follows from the rela- 
tionship between the eigenvector pi and the Uj being pi = {ai — 02)/ V^ in our chosen phase convention. 



III. (UN)RELATING FLAVOUR SYMMETRIES TO THE Upmns FORM 



This formalism can be used to systematically explore symmetries that are consistent with TB form. 
Consider an exact {TB) form of Upmns at leading order in the FSE. This is dictated by 613 = and 
^23 maximal. We restrict the Yukawa vectors in this discussion to real values for simplicity. Including an 
unfixed S12 this gives (in a particular phase convention) 



/ 



Utb 



C12 S12 



\ 



^12 C12 
V2 \f2 



1 
V2 



(29) 



\ _£i2 £12 J_ I 

\ V2 V2 V2 / 
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Assuming Upmns = Utb 
1 



1 



C12 ,^ , ^ . 

P2 = ^(CTl +0-2J + Sl2 0"3, 

P3 = -^(i?l + (T2) +C12(?3, O's = S12P2 + C12/O3. (30) 

Consider the "flavour" basis as an illustrative example of the FSE, where ffi = (0,0,1), (T2 



Tl = -j={pi+ Cl2 P2 - Sl2 Ps), 

0-2 = -y^{-Pi + C12 /02 - S12 ps), 

O's = S12P2 + C12/O3. 



(0,1,0), (T3 = (1,0,0). In this case, Z//(e, L) = diag{l, 1,1) and Upmns =U{v,L). The flavour ba- 
sis is appealing in that it allows the solution for the perturbations to the pi eigenvectors and the Upmns 
matrix to proceed simply through solving 



-XX + - — yy 



'^fJ'a 



2/Ufc' 



U^p diag(0, 5mb, nia + 6ma)ul^, 



(31) 



/ 



-. -.T , pa -. -cp 

XX H yy 



{rUa + 5ma) Pa 







\ 



+ 



\/2 5mhPa 



I \/2s^2 ■512 Cl2 SI2C12 » 



S12 C12 



■^12 



"^12 



s/2 s/2 



I 



\^i2ci2 ^ ^ ; 



1 -1 

yo -1 1 

Trivially one finds 

(0, -1, 1) ^J{ma + bma) Pa/v, if = f\/2si2,ci2,ci2J y^6mbPb/v. (32) 



X 



It follows that in the flavour basis x ■ y = so 5ma = and ((5p^, 5ph, 6pc) = (0, 0, 0). Now consider 
including a third neutrino eigenvalue, retaining the required perturbation. For a solution, Z2 = 2:3 is required, 
and consequently 62171^ = as x • i* = 0. As a result the leading eigenvector pa is unperturbed by integrating 
out both Nh, Nc in the flavour basis, one finds 

^/2p^^S2rnc 



z^ = ±- 



cf,+ 



(527716 2 1 / (^27716 



C-19 + Si 



1 /^2777fe 



C10 + Si 



(33) 



^12 ^ ^2777, ^12' ^ Y 5^^^ -12 - ^12, ^ Y 52777^ ^^^ ^ ^12^ ' 

A /i -H- r symmetry implemented on yand zis consistent with TB form of the PMNS matrix, as expected. 



A. Perturbative breaking of TB form 

Now consider the breaking of TB form. The value of ^13 measured by the DAYA-Bay collaboration is 
in agreement with the global fit values given in Ref. 131, as such, we use the fit results to find the measured 
pattern of deviations in TB form. It is instructive to construct the following ratios of experimental values 
characterizing the deviations of TB form in each mixing angle. Using the small angle approximation 



tan^(6ei2) tan^(6e23) tan^(6ei2) 

■ 2) J =0.02 ±0.32 ■ 2) J =0.01 ±0.10 , ^' 

sm"^ ((56*13) sin"^(<56'i3) tan^((56'23) 



2.0 ±36. 



(34) 
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Here we have used the one sigma new reactor flux values of Ref. lO and taken a it symmetric one sigma 
error. Various breaking of TB form can be studied using the FSE and compared to these results. Consider 
the case that Ni, retains a fi ■^ t flavour symmetry in its couplings to the charged leptons, but Nc breaks 
such a symmetry so that deviations in TB form are expected. Fix z^ = z^ + (0, Ai, A2) with Ai 7^ A2 
and treat this breaking as a perturbation using the FSE. We use Am\^ ~ Am\(j assuming a normal 
hierarchy. At leading order in the FSE the breaking of TB has the pattern 

tan^(J6'23) _ 2 cos^ Op^^ + cos^ 9p^^ 
sin^ (SOis) cos2 Op^^ + 2 cos^ Op^^ ' 

KtI^ = »• 0« 

tan^(o6'23) 
As the range of the predicted value of tan^ (6623) / sin^ {6613) is given by [0.5, 2] at leading order in the 
FSE, this pattern of flavour breaking does not reproduce the data. In this way, particular mechanisms of the 
breaking of TB form can be falsified. One can perform the exercise of assuming TB form is broken in the 
flavour basis by Nb so that y^ = 2/^ + (0, Ai , A2). Using the expansion one finds the pattern of deviations 
are distinct. These breakings of TB form are also correlated with mass perturbations in each case which 
are trivial to determine using this formalism. 

We emphasize however that the relationships between the eigenvectors determine the form of the PMNS 
matrix in general. This is easy to demonstrate in more detail. Consider retaining a /i -H- r symmetry 
imposed on y and zbut deviating from the flavour basis, choosing x, y, z as above but unfixing cjj. At 
leading order in the expansion of Upmns one can solve for the condition that TB form is still obtained 
for general ai. One finds that only the flavour basis for ffi gives a valid solution at leading order in the 
expansion. This makes clear that ^ ■(r^ t symmetry imposed on the Lagrangian alone is not related to TB 
form in a basis independent manner. 

As a further example that fi ■(r^ t symmetry is also not unique or of particular physical significance 
in allowing TB form (with an appropriate choice on the Ui eigenvectors), consider the following pro- 
cedure. Choose a {e,fj,,T) symmetry on the first interaction eigenvector, ie x^ = (l,l,l)/-v/3, and 
(1, 1,0) /\/2 as a simple interaction eigenvector for Ni, leading to an orthonormal eigenbasis at 



dP 



y 



leading order, finding 



U{u,L) 



/ 


1 

V2 


1 1 \ 




1 


1 1 




V2 


V6 V3 


V 





~V 3 73 / 



(36) 
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Solving directly for U{v, L) so that at leading order TB form is obtained, one finds 

^1(^/2 + 2^3) i (^/2 - V3 - ^/6) 1(^2-^/3-^/6)^ 

U{e,L)= 1(^-2^3) i(V2 + V3-^/6) 1(^2 + ^3 + ^/6) • (37) 

1 _ V2 2+x/3 _ V2 I J_ 

\ 3 3^2 3 "^ V6 / 

This procedure can be used for any flavour symmetry chosen to fix x, y in the FSE for the Ni. 

We also note that the impact of sterile neutrinos weakly coupled to the SM on neutrino phenomenology 
can be systematically studied with this approach. For example, one can relate any measured value of a 
deviation from TB form to the particular- Yukawa coupling vector of a single sterile neutrino, which can 
be shown to accommodate the value of ^13 reported by the Daya-Bay collaboration while the three right 
handed neutrinos partners of the SM fields give an exact TB form of the Upmns matrix. 

IV. CONCLUSIONS 

Flavour symmetries that are related to the structure of the Upmns matrix only in a particular basis 
choice ofU{e, L) can lead to suspect physical conclusions. As an alternative to basis dependent symmetry 
studies, we have developed a perturbative expansion relating the measured masses of the neutrinos to the 
form of the PMNS matrix. This expansion offers a promising framework for broadly understanding the 
implications of the systematically improving experimental neutrino data, particularly in a normal hierarchy. 

We have illustrated our approach in an example where the flavour basis was chosen, for the sake of 
familiarity, and then shown how the expansion can control the predictions of TB form being broken. How- 
ever, the approach we outline can accommodate any basis choice. Indeed, it is the relationships between 
the eigenvectors that dictate the form of the Upmns matrix in a U{e, L) basis independent manner. This 
formalism can be employed in model building to attempt to determine a compelling origin of the eigenvec- 
tor relationship that corresponds to TB form at leading order in the FSE. Further, as the breaking of TB 
form is now experimentally established due to the discovery of a non zero ^13 in Ref. [7], we expect the 
FSE to be of some phenomenological utility in falsifying mechanisms of the breaking of the TB form of 
the UpMNS matrix, as the pattern of this breaking is further resolved experimentally in the years ahead. 
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